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Effect of disorder on superconductivity 
in the boson-fermion model 
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We study how a randomness of either boson or fermion site energies affects the superconducting 
phase of the boson fermion model. We find that, contrary to what is expected for s-wave super- 
conductors, the non-magnetic disorder is detrimental to the s-wave superconductivity. However, 
depending in which subsystem the disorder is located, we can observe different channels being af- 
fected. Weak disorder of the fermion subsystem is responsible mainly for renormalization of the 
single particle density of states while disorder in the boson subsystem directly leads to fluctuation 
of the strength of the effective pairing between fermions. 



PACS numbers: 74.20.-z, 74.20.Mn, 74.25.Bt, 71.10. 
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I. INTRODUCTION 

The boson fermion (BF) model is an example of a 
microscopic theory of nonconventional superconductiv- 
ity. It describes a mixture of itinerant electrons or holes 
(fermions) which interact via charge exchange with a sys- 
tem of immobile local pairs (hard-core bosons). Due to 
this interactions, bosons acquire finite mass and under 
proper circumstances might undergo Bose condensation 
transition while fermions simultaneously start to form a 
broken symmetry superconducting phase. 

For the first time this model has been introduced ad 
hoc almost two decades ago to describe the electron 
system coupled to the lattice vibrations in a crossover 
regime, between the adiabatic and antiadiabatic limits. 
Later it has been formally derived from the Hamiltonian 
of wide band electrons hybridized to the strongly corre- 
lated narrow band electron system || . Very recently || 
the same effective BF model has been derived purely from 
the two dimensional Hubbard model in the strong inter- 
action limit using the contractor renormalization method 
of Morningstar and Weinstein [Q] . 

Some authors have proposed it as a possible scenario 
for description of high temperature superconductivity 
(HTSC). The unconventional way of inducing the super- 
conducting phase in the BF model has been indepen- 
dently investigated in a number of papers [f^|il"|] . More- 
over, this model reveals also several unusual properties 
of the normal phase (for T > T c ) with an appearance of 
the pseudogap being the most transparent amongst them 
Jl2| pd[ . Apart of eventual relevance of this model to 
HTSC there are attempts to apply the same type of pic- 
ture for a description of the magnetically trapped atoms 
of alkali metals 

The important question which we want to address in 
this paper is: what is an influence of disorder on super- 
conductivity of the BF model ? The conventional s-wave 
symmetry BCS-type superconductors are known to be 



rather weakly affected by paramagnetic impurities - 
the fact which is known as "Anderson theorem". Non- 
magnetic impurities have remarkable detrimental effect 
on superconductors with the anisotropic order parame- 
ters. Magnetic impurities lead to pair-breaking effects 
which result in a strong reduction of T c even in s-wave 
superconductors. Studying the effect of impurities on 
the superconductors has always been an established tool 
for investigation of the internal structure of the Cooper 
pairs. 

Due to the nonconventional pairing mechanism (i.e. ex- 
change of the hard-core bosons between fermion pairs) it 
is of a fundamental importance to see how the nonmag- 
netic impurities (disorder) affect the isotropic supercon- 
ducting phase of the BF model. Previously, such a study 
has been carried out by Robaszkiewicz and Pawlowski 
|l7f who considered disorder only in the boson subsystem. 
Using a method of configurational averaging for the free 
energy, authors have shown a strong detrimental effect 
of disorder on superconductivity. Apart of a reduction 
of the transition temperature T c they have also reported 
a remarkable change of a relative ratio A(T = 0)/fcsTc. 
In this paper we analyze the effect of disorder present 
in both: fermion and boson subsystems using a different 
method of the coherent potential approximation. 



II. THE MODEL AND APPROACH 

A. Hamiltonian of the disordered BF model 

We consider the following Hamiltonian of the disor- 
dered BF model 

H BF = ^ tijclrCjv + ^ ( £ * ~ ^) 4* c «? 
i,j.(T i 

+ (A B + Ek - 2fx) b\b t 



1 



(i) 



C. Fermion part 



We use the standard notations for annihilation (creation) 
operators of fermion Cj )Cr (cj cr ) with spin cr and of the hard 

core boson fcj (&J) at site i. Fermions interact with bosons 
via the charge exchange interaction v which is assumed 
to be local. There are two ways in which disorder enters 
into the consideration. Either (a) fermions are affected 
by it and this is expressed by the random site energies 
£i, or (b) hard core bosons via their random site energies 
Ei, 

To proceed, we apply first the mean field decoupling 
for the boson fermion interaction 



b l*i*i - (h)* c 4 c *t + b t (cii^r) 



(2) 



which is justified until v is small enough in compari- 
son to the kinetic energy of fermions. After decoupling 
(||) we have to deal with the effective Hamiltonian com- 
posed of the separate fermion and boson contributions 
H ~ H F + H B f|§ 



H 



E i -2 f i) b\bi + Xi 



(3) 



(4) 



where x^ = v (c^c^) and pi = v(bi). The site depen- 
dence of pi and Xi indicates the disorder induced ampli- 
tude fluctuations of the order parameters. 



B. Boson part 

For a given configuration of disorder we can exactly 
find the eigenvectors and eigenvalues corresponding to 
the lattice site i using a suitable unitary transformation. 
Statistical expectation values of the number operator b\bi 
and the parameter pi are given by 



bib, 



l__A B+Ei -2, p, 
2 4 7l \k B T)' 



vxi , I 7, 
— - tanh — — 

2 7 * \k B T 



(5) 
(6) 



where 7$ = j\/(Ab + E'i — 2^) 2 + 4\xi\ 2 and /c^ is the 
Boltzmann constant. Note, that the site dependent 
fermion order parameter Xi enters the expression for the 
boson number operator (^|) and the parameter pi 
Disorder of any subsystem is thus automatically trans- 
fered onto the other one. 



Analysis of the fermion part (y) is more cumber- 
some. To study it we use the Nambu representation 
\&t = (cL, Cii), #i = (^lY and define the matrix Green's 

function G(i,j;uj) = ^/^j;^t\^ . Equation of motion 
for this function reads 



E 

i 



(w - si + p)5 t .i - Ui ~P*hi 

-Pi8 l i (uj + ei - p)8 a + Ui 

x G(l,j;uj) = 1%. 



(7) 



Using the matrix Green's function G°(i, j; uj) of a clean 
system 



[G°(k; 



£k + P 

— p uj - 



-P" 



(8) 



where p = vj^ (pi) is a global order parameter (which 
plays a role of the effective gap in the superconducting 
fermion subsystem) , and defining the single site impurity 
potential Vi as 



Vi = 



-Pi 



-pi 

-si 



(9) 



one can write down the following Dyson equation for the 
Green's function G(i,j;uj) 

G(i, j; uj) = G°(i, j; uj) + J2 G°(i> h «)V,G(i,i; uj). 

i 

(10) 

This Green's function depends on the specific disorder 
configuration. In order to pass through one usually av- 
erages it over the all possible configurations. 

For carrying out the configurational averaging we use a 
method of the Coherent Potential Approximation (CPA). 
The main idea of CPA is to replace the random potential 
Vi by some uniform coherent potential S(cj). Formally, 
the Green's function which satisfies ( fio| ) with V/ replaced 
by E(w) is then given (in the momentum coordinates) by 



G 



CPA 



(k;u;) 



[G°(k; W )] X -E(u;). 



(11) 



Configuration at site i is defined by values of the ran- 
dom energies £j, E{ - we shall symbolically denote it by 
a = {ei, Ei}. Any of possible configurations a can occur 
with some probability P({ei, El }) = c (a \ and of course 
these probabilities are normalized ^2 a c^ = 1. 

A particle propagating through the medium character- 
ized by the coherent potential £(u;) is thus, at site i, scat- 
tered with probability by the potential — S(w). 
For a chosen configuration a of the site i the conditionally 
averaged local Greens function is given by 



G 



CPA 



V 



(a) 



E(w) 
(12) 



This Green's function G^ a \i,i;uj) describes the system 
in which all sites, except one indicated by i, are described 
by the coherent potential S(w). In CPA one requires 
that, the average of the local Green's function is the same 
as the Green's function of the averaged system. This 
CPA condition is identical with the following equation 



(13) 



Equations (JTT|-jT^) have to be solved selfconsistently 
to yield the coherent potential £(w). Physical quantities 

such as fermion concentration n F = jj J2i a i^ c \<j c i<^j an< ^ 
the superconducting order parameter x = ■h ^ x % are to 
be calculated from 



1 lm[Gg PA (t, 1 ;cj + t + )} (15) 



e /3w 
dbJ 



2 
ttJV 
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where j3 = l/k B T. 

In the following section we discuss the changes of the 
superconducting transition temperature T c caused by dis- 
order. 



III. DISORDER IN FERMION SUBSYSTEM 

It is instructive to investigate the disorder separately 
for fermion and boson subsystems. Let us start with 
fermion disorder £j. We set = for all lattice sites. 
For the random fermion energies we choose Si = Sq with 
probability c and = with probability 1 — c. It is a 
bimodal type disorder 



P ({£,}) = cSie, - e ) + (1 - c) 5{e l - 0). 



(16) 



Here we shall be mainly interested in the supercon- 
ducting transition temperature T c . In this limit ]l9| the 
diagonal disorder affects mainly a diagonal part of the 
matrix Green's function G. In fact, even for no disorder 
acting directly in bosonic subsystem the boson order pa- 
rameter p in equation (j^) does depend on the site index 
via fermion order parameter Xi. However, we expect this 
induced disorder to be weak and neglect it. This allows us 
to show how disorder in fermionic subsystem only, affects 
T c . 

The off-diagonal elements of the coherent potential 
vanish. Due to the general symmetry £22 (iw) = 
— £n (— iuj) fl9|| we can simplify the self-energy matrix 
to 



£(«j) = 



Sii(iw) 






-£ n (-Mj) 



(17) 



£11(0;) can be found from the CPA equation (13) which, 
for a normal phase, takes a well known form [18| 

(18) 



[Euiuj)}- 1 + F(cj) [Eii(w)-eo] 1 + F(cj) 

with F(uj) = ^Ek G n PA ( k ' w )- Equation Q should 
be solved subject to a given dispersion relation and 
parameters c, Eq. 

Finally having calculated £n(w), we can find n F and 
x ( |T4| , |l5| ) as well as n B , p (Spj). In particular, the critical 
temperature T c = {ksPc) is given via 

2 tanh \p e (A B - 2/x) /2] 



1 



2(i 



E 



dwi 



x A(k, ^(k, . 2 ) tanh + tmh [M2] (19) 

2 (LOl + LU 2 ) 

where A(k.,w) = (-l/7r)Im | G^ PA (k, uj + i0+)\ de- 
notes the spectral function of the normal phase. 

We choose for our study a case of weak boson fermion 
interaction v = 0.1 (in units of the initial fermion band- 
width) and total concentration of charge carriers n to t = 
2ns + np = 1 . Figure ([j]) shows how T c of a clean sys- 
tem depends on position of the boson level. There are 
three distinguishable regimes J^,|| of relative occupancy 
by bosons and fermions. Superconducting correlations 
are of course most visible when chemical potential is close 
to Ag/2. We choose the value A^/2 = —0.3 to be close 
to optimal value of transition temperature and to have 
comparable amount of fermions and bosons. For compu- 
tations we use the 2D square lattice dispersion - the van 
Hove singularity is safely distant from the Fermi energy 
for the above parameters. 
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FIG. 1. Variation of T c with respect to boson energy As 
for a clean system with ntot = 1- Bottom panel illustrates the 
concentrations of fermions (n F ) and bosons (n B ) at T = T c . 
Note the three distinct regimes of: predominantly local pairs 
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n to t, coexisting pairs and fermions n 



2n B , and 



predominantly fermions n ~ ntot (so called BCS limit). 



In figure g we plot the transition temperature T c , cal- 
culated from equation (|l9|) against concentration c for 
several values of £o. With an increase of concentration 
c of scattering centers we notice a gradual reduction of 
the critical temperature. This tendency can be under- 
stood by looking at the behavior of the fermion density 
of states at the Fermi energy q(sf)- Disorder is responsi- 
ble for renormalization of the low energy sector and these 
low energy states are involved in forming the supercon- 
ducting type correlations. As shown in the bottom panel 
there is additional effect coming from the rearrangement 
of occupations n F and n B . With an increasing concentra- 
tion c the fermion band is shifted toward higher energies 
and the system is then mainly occupied by bosons (so 
called, the local pair LP limit). 
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FIG. 2. Transition temperature T c as a function of concen- 
tration c of scattering centers with various positive values of eo 
(top panel). Density of states at the Fermi energy g{eF) (mid- 
dle panel) and relative occupations by bosons and fermions 
(bottom panel) for eo = 0.5 at T = T c . 



For negative values of £o the disorder shows stronger 
influence on T c . On one hand we have again a direct 
effect of the renormalized density of states (see g(ep) in 
the middle panel of figure 3). On the other hand, with 
an increase of c for any negative value of Eq the fermion 
band and the position of the chemical potential drift to- 
wards lower energies. As its consequence the number of 
fermions increases and the number of bosons decreases. 
Effectively we thus approach the BCS limit where transi- 
tion temperature diminishes very fast if Ab /2 goes above 
/i (check for example the curves for £o = —0.4 and —0.5). 
The strong disorder in fermion subsystem makes the pair- 
ing mechanism almost ineffective at all. 



0.015 

T. 



c 
0.01 



0.005 



1 

0.5 



0.5 



goZ-0-1 



Nr0.3 



\ V0.4 
-0.5 \ \ 



9(e F ) 




0.2 0.4 0.6 0.8 q 1 

FIG. 3. The same as in figure (^) except that for nega- 
tive values of eo (top panel). The middle and bottom panels 
correspond to eo = —0.5. 

In figure | we plot T c versus (positive) £o for several 
concentrations c. Again, T c roughly follows variation of 
the density of states 3(cf) which is shown in the bottom 
panel. As discussed above for large values of concentra- 
tion c and large positive £o the system is mainly filled 
by bosons (the LP limit) so there is some finite T c even 
when fi — As/2 is far below the fermion band, this is an 
artifact of the mean field approximation |2j,p| . 
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FIG. 4. Transition temperature T c as a function of the 
energy eo of the scattering centers whose concentration is c 
(top). Density of states <7(cf) for each of the concentrations 
c (bottom panel) For c = 1 and for eo > 0.2 Fermi energy 
goes below the fermions band, system is then strictly in the 
LP limit of the BF model. 

Behavior of T c with respect to negative values of £o 
can be easily deduced from the figure 3 so we skip this 
illustration. 

In summary we notice that change of the transition 



A 



temperature T c caused by weak disorder in fermion sys- 
tem is controlled mainly by modification of the low lying 
energy states. This is in accord with the Anderson the- 
orem for spin singlet s-wave superconductors. However, 
additional influence comes from redistribution of particle 
spectrum and their relative occupancy and such effects 
arc dominant for large values of impurity concentration c 
and for their large scattering strength |eo|- In this limit 
the boson - fermion exchange becomes ineffective. 



IV. DISORDER IN BOSON SUBSYSTEM 

Now we turn attention to a case when boson energies 
are random Ei ^ and, for simplicity, assume no fermion 
disorder i.e. Si = for all the lattice sites. The scattering 
potential (S) reduces then to 



V, 



o -fi(cU 



-// (c/|Q T ) 



with 



and 



fl=v 



2 tanh [flji] 



2 7 i 



'Ab + Ei 



(20) 



(21) 



(22) 



It means that the fluctuating boson energy level Ei in- 
duces fluctuations of the pairing strengths /; in the 
fermion subsystem. To some extent, this situation re- 
minds the negative U Hubbard model |2(]] for which the 
random local attraction Ui < leads to the following 
scattering matrix 



Hub) 



Ui(ni)/2 Uiic]^) 
Ui{c Ll c n ) Ui(ni)/2 



(23) 



We see that in our case the role of a random pairing 
potential U\ is played by —j\ given in equation (pi]). 

There are two extreme limits, as far as the effectiveness 
of the random boson energy Ab + Ei is concerned 

• for small (on the scale of fermion-boson interac- 
tion v) fluctuations of E h effect of the disorder 
becomes negligible unless the chemical potential is 
pinned to the boson level fi = As/2, when the 
amplitude of the pairing potential is controlled by 
fi ~ w 2 tanh \(3x{\ /2xi and is usually uniform except 
at very low temperatures (3 — > oo when fi ~ v 2 jxi 

• for large fluctuations of Ei one obtains fi ~ 



w 2 tanh 



f (A B + Ei- 2/i) I {A B + E t - 2p). 



To analyze effects of the disorder in boson sub- 
system we use a two pole distribution P({Ei}) = 
\ [S(Ei - E ) + S(Ei + E )}. The boson energy is A B ± 
-Eo with an equal probability 0.5. Figure (||) shows criti- 
cal temperature T c , calculated from equation ([l5|), as a 
function of energy Eg by which the boson energy is split. 
Strong dependence of T c on disorder is a combined effect 
of the density of states, the fluctuating interactions and 
the changes in concentration of carriers. 
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FIG. 5. Transition temperature T c (top), the averaged 
pairing potential < /; >= J2{E t } P i{ E i))fi (middle), to- 
gether with the occupation of fermions rip and bosons ns 
at T = T c (bottom). 
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FIG. 6. Normalized critical temperature T c /T c (Eo = 0) 
and the normalized pairing potential </;>/< /; (-Bo = 0) > 
versus energy Eq. T^ bcs ^ shows the BCS-like relation be- 
tween critical temperature and pairing potential. Left panel 
refers to ntot = 1, Ab = —0.6 discussed above and the right 
panel corresponds to the symmetric case of the BF model 
A_b = 0, ntot = 2 (with the half filled boson and fermion 
subsystems) . 

To estimate what influence comes only from the renor- 
malization of the effective pairing we plot in figure | the 
normalized transition temperature denoted by T c and the 
normalized averaged < fi > for the parameters given 
above (left panel), and for a fully symmetric case of 
the BF model (right panel). The transition temperature 
T c is the BCS-type estimate of the effect of changes 
in the effective pairing due to disorder 

T C (E ) cx exp ( ~\ \ . (24) 

A general trend observed in figure || is that the average 
effective interaction < Ji(Eq) > decreases with increasing 
disorder, even though the bar fermion-boson interaction 
v remains constant. This decreasing pairing interaction 
is the only factor responsible for a behavior of T c versus 
Eq in the right panel. We notice absence of the BCS- 
like exponential scaling which is due to unconventional 
pairing in the BF model. 

In the left panel, corresponding to the above studied 
case Ab — —0.6, n tot — 1, we notice a larger discrepancy 
between the pairing amplitude and T c . With an increase 
of Eq the transition temperature is much strongly re- 
duced than in the symmetric case. This effect has to be 
assigned to redistributions of particle occupancies. At 
large values of Eq we have practically only hard core bo- 
son particles in the system, and they cannot induce su- 
perconductivity among fermions whose fraction becomes 
very small. 

In the previous study jn| authors have used the same 
bimodal distribution of random boson energies. They 
have found a strong reduction of T c near Eq ~ 2v which 
agrees well with our data shown in figure 6. Moreover, 
the authors have reported that disorder affects the ra- 
tio A SC (T = 0)/fcsT c which changes from 4.2 (for a clean 
system to the standard BCS result 3.52 at large Eq. 
Simple explanation of this effect can be offered. The bo- 
son energy (which is split by 2Eq) is for sufficiently large 
Eq partly in the LP limit (for Ei — —Eq) and partly in 
the BCS limit (if Ei = +Eq). The second one contributes 
with the standard BCS value if |2?o| is large enough (see 
e.g. Fig. 9 in Ref. f§). 

V. CONCLUSION 

The randomness of the site energies of both, fermions 
and bosons, has a strong effect on superconducting phase 



of the BF model. Weak disorder in the fermion subsys- 
tem affects the superconducting transition temperature 
mainly via rescaling the low energy states which are in- 
volved in the the formation of the Cooper pairs. There- 
fore T c roughly follows the density of states at the Fermi 
energy. For sufficiently large disorder Eq there appears 
some critical concentration c at which T c may eventually 
drop to zero. 

Disorder in the boson subsystem has a much more fine 
influence on superconductivity. Randomness of boson 
energies is transformed directly into randomness of the 
pairing strength. Effectively physics of the disordered 
BF model becomes similar to that of the random neg- 
ative U Hubbard model [^0|. Even the relatively small 
fluctuations of the boson energies show up their strong 
detrimental effects on superconductivity. 

In a simple minded picture one can envision this sit- 
uation as a random change between various regimes of 
superconductivity. Depending on a value of E\ the bo- 
son energy Ab + Ei can be either far below the Fermi 
energy (the LP limit), or far above the Fermi energy 
(the BCS limit). Each of such random configurations 
contributes with a different strength of superconducting 
correlations. On average, the superconducting transition 
temperature T c strongly diminishes and practically dis- 
appears if the amplitude of the randomly fluctuating bo- 
son energies \Eq\ is large enough. 

In summary, our calculations show that disorder 
strongly affects the s-wave superconducting phase of the 
BF model. This apparent contradiction with Anderson 
theorem can be understood because of a change of the ef- 
fective pairing interaction induced by disorder, and this 
effect is contrary to the Anderson's main assumption |l6| . 

To compare our results with experimental data on high 
temperature superconductors one has to consider the d- 
wave superconducting order parameter. This type of a 
symmetry arises in a natural way according to the re- 
cent derivation of the BF model . Effect of disorder on 
such anisotropic superconducting phase of the BF model 
is outside the scope of the present paper and will be dis- 
cussed elsewhere. 
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